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Abstract 

Structure of exact solutions modelling confinement is discussed for SU(3)-Yang-Mills 
equations and uniqueness of such confining solutions is proved in a certain sense. 
Relationship of the obtained results to QCD and the confinement mechanism is 
considered. Incidentally the Wilson confinement criterion for the found solutions is 
verified and also numerical estimates for strength of magnetic colour field responsible 
for linear confinement are adduced in the ground state of charmonium. 
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1 Introduction 

As was remarked in Ref. [1], there exists a natural way of building meson 
spectroscopy and relativistic models of mesons which might be based on the 
exact solutions of the SU (3)- Yang- Mills equations modelling quark confine- 
ment, the so-called confining solutions, and also on the corresponding modulo 
square integrable solutions of the Dirac equation in those confining SU(3)- 
fieids. The given approach is the direct consequence of the relativistic QCD 
(quantum chromodynamics) Lagrangian since the mentioned Yang-Mills and 
Dirac equations are derived just from the latter one. In Ref. [1] both the types 
of solutions were obtained and then in Ref. [2] they were successfully applied 
to the description of the charmonium spectrum. In its turn, the mentioned 
description pointed out the linear confinement to be (classically) governed by 
the magnetic colour field linear in r, the distance between quarks. 
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The results of Refs. [1,2] suggested the following mechanism of confinement to 
occur within the framework of QCD (at any rate, for mesons and quarkonia). 
The gluon exchange between quarks is realized in such a way that at large dis- 
tances it leads to the confining SU(3)-field which may be considered classically 
(the gluon concentration becomes huge and gluons form the boson condensate 
- a classical field) and is a nonperturbative solution of the SU (3)- Yang-Mills 
equations. Under the circumstances mesons are the relativistic bound states de- 
scribed by the corresponding wave functions - nonperturbative modulo square 
integrable solutions of the Dirac equation in this confining SU(3)-field. For 
each meson there exists its own set of real constants (for more details see Sec- 
tion 2) aj, Aj,bj, Bj parametrizing the mentioned confining gluon field (the 
gluon condensate) and the corresponding wave functions while the latter ones 
also depend on /i , the reduced mass of the current masses of quarks forming 
meson. It is clear that constants aj, Aj,bj, Bj, /i should be extracted from 
experimental data. Further application of the approach to quarkonia (char- 
monium and bottomonium) in Refs. [3,4] confirmed the mentioned picture of 
linear confinement and, in particular, gave possibility to estimate the above 
gluon concentrations. 

The aim of the present Letter is to specify a number of features in the above 
confinement scenario, in particular, to show that the confining solutions found 
in Ref. [1] and used in Refs. [2-4] are in essence the unique ones. 

Further we shall deal with the metric of the flat Minkowski spacetime M that 
we write down [using the ordinary set of local spherical (r, i?, ip) or rectangular 
(Cartesian) (x,y, z) coordinates for spatial part in the forms 

ds 2 = g^dx" ® dx v = dt 2 - dr 2 - r 2 (d$ 2 + sin 2 $d^ 2 ) , (1) 



or 

ds 2 = g^dxP <g> dx v = dt 2 - dx 2 - dy 2 - dz 2 , (2) 

so the components g^ u take different values depending on the choice of co- 
ordinates. Besides we have 5 = \det(g fiu )\ = (r 2 sini9) 2 in spherical coordi- 
nates and the exterior differential d = d t dt + d x dx + d y dy + d z dz or d = 
d t dt + d r dr + d^dd + d v df in the corresponding coordinates and we denote 
3-dimensional vectors by bold font. 

Throughout the paper we employ the system of units with h = c = 1, unless 
explicitly stated otherwise. 
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2 Uniqueness of the confining solutions 



If A = A^dx^ = A^Xadx^ is a SU(3)-connection in the (trivial) three-dimensional 
bundle £ over the Minkowski spacetime, where A c are the known Gell-Mann 
matrices, then we are interested in the confining solutions A of the SU(3)- 
Yang-Mills equations 



while the curvature matrix (field strentgh) for the ^-bundle is F = dA + gA A 
A = F^ u \ a dx^ A dx v and * means the Hodge star operator conforming to 
metric (1), g is a gauge coupling constant. 

The confining solutions were defined in Ref. [1] as the spherically symmetric 
solutions of the Yang-Mills equations (3) containing only the components of 
the SU(3)-field which are Coulomb-like or linear in r. Additionally we shall 
impose the Lorentz condition on the sought solutions. The latter condition is 
necessary for quantizing the gauge fields consistently within the framework 
of perturbation theory (see, e. g. Ref. [5]), so we should impose the given 
condition that can be written in the form div(A) = 0, where the divergence 
of the Lie algebra valued 1-form A = A^dx^ = A®\ a dx^ is defined by the 
relation (see, e. g. Refs. [6]) 



It should be noted the following. If writing down the Yang-Mills equations (3) 
in components then we shall be drowned in a sea of indices which will strongly 
hamper searching for one or another ansatz and make it to be practically 
immense. Using the Hodge star operator as well as the rules of external calculus 
makes the problem to be quite foreseeable and quickly leads to the aim. Let 
us remind, therefore, the properties of Hodge star operator we shall need (for 
more details see Refs. [6]). 

2. 1 Hodge star operator * on Minkowski spacetime in spherical coordinates 

Let M is a smooth manifold of dimension n so we denote an algebra of smooth 
functions on M as F(M). In a standard way the spaces of smooth differential 
p-forms A P (M) (0 < p < n) are defined over M as modules over F(M). If a 
(pseudo)riemannian metric G = ds 2 = g^dx^ <g> dx v is given on M in local 
coordinates x = (x M ) then G can naturally be continued on spaces A P (M) by 
relation 



d* F = g(*F A A — A A *F) , 



(3) 



div(A) = JLd^Vs^) . 



(4) 



G(a,0) = det{G(a h fa)} 



(5) 
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for a — a-y Aa 2 ... Aa p , j3 — fli A(3 2 ... Aj3 p , where for 1-forms = a^dx 11 , f3j = 
Pi j) dx u we have G{ai ) f3j) = g^a^}^ with the Cartan's wedge (external) 
product A. Under the circumstances the Hodge star operator *: A P (M) — > 
A n ~ p (M) is defined for any a e A P (M) by 

a A (*a) = G(a, a)u g (6) 



with the volume n-form oj g = y | det^^cte 1 A...dx n . For example, for 2-forms 
F = F^dx* 1 A fir 1 ' we have 

F A*F = (g^g^ - g^ p g ua )F^F aP ^~5 dx l Adx 2 ---A dx 11 , n < v, a < (3 (7) 

with 5 = | det(g M j,)|. If s is the number of (-1) in a canonical presentation of 
quadratic form G then two most important properties of * for us are 

* 2 = (_1)P("-P)+* , ( 8 ) 

+ /2«2) = + h(*oi 2 ) (9) 

for any /i, / 2 G F(M), a±, a 2 € A P (M), i. e., * is a F(M)-linear operator. Due 
to (9) for description of *-action in local coordinates it is enough to specify 
*-action on the basis elements of A P (M), i. e. on the forms dx n Adx 12 A...Adx lp 
with i\ < i 2 < ... < i p whose number is equal to C p = ( n _^) !p! ■ 

The most important case of M in the given paper is the Minkowski spacetime 
with local coordinates t, r, ip, where r, ip stand for spherical coordinates on 
spatial part of M . The metric is given by (1) and we shall obtain the *-action 
on the basis differential forms according to (6) 

*dt = r 2 sin §dr A dd A dip, *dr = r 2 sin fldt A d-d A dip, 

*dd = — r sin ddt Adr A dip, *dip = rdt A dr A dd, 
*(dt A dr) = -r 2 sin fidd A dip , *(dt A dfi) = sin fidr A dip , 

*(dt A dip) = — r -dr A dfi , *(dr A dti) = sin fidt A dip , 

sinw 

*(dr A dip) = — t dt /\ dfi /\ dip) = — -dt A dr , 

sin v r 2 sin v 

*(dt Adr A d§) = -dip , *(dt A dr A dip) = -d-d, 

r r sin v 

*(dt Add A dip) = - 1 - dr, *(dr A dd A dip) = - 1 dt, (10) 
r 2 sin v H sin i? 

so that on 2-forms * 2 = — 1, as should be in accordance with (8). 

At last it should be noted that all the above is easily over linearity continued on 
the matrix- valued differential forms (see, e. g., Ref. [7]), i. e., on the arbitrary 
linear combinations of forms a ili2 ...i p dx n A dx' 12 A ... A dx %v , where coefficients 
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Oiii 2 ...ip belong to some space of matrices V, for example, a SU(3)-Lie alge- 
bra. But now the Cartan's wedge (external) product A should be understood 
as product of matrices with elements consisting of usual (scalar) differential 
forms. In the SU(3)-case, if T a are matrices of generators of the SU(3)-Lie 
algebra in 3-dimensional representation, we continue the above scalar product 
G on the SU(3)-Lie algebra valued 1-forms A = A a ^T a dx p and B = B b v T b dx v 
by the relation 

G(A,B)=g^A;B b u Tr(T a T b ), (11) 

where Tr signifies the trace of a matrix, and, on linearity with the help of (5), 
G can be continued over any SU(3)-Lie algebra valued forms. 



2.2 The confining solutions 



Let us for definiteness put T a = \ a . Under this situation we can take the 
general ansatz of form 

A = rTdt + A r dr + A 6 d& + T v My, (12) 

where A$ = r p T and matrices f = a a \ a , A = /3 a X a , T = 7 a A a , A r = 
f a (r)\ a with arbitrary real constants a a ) f3 a ) ^j a and arbitrary real functions 
f a (r). It could seem that there is a more general ansatz in the form A = 
r Pa a a \ a dt + f a (r)X a dr + r par ) a \ a d'd + r Va f3 a \ a dip but somewhat more compli- 
cated considerations than the ones below show that all the same we should 
have jj, a = fJ>,v a = u, p a = p for any a so we at once consider this condition 
to be fulfilled to avoid unnecessary complications. Then for the above ansatz 
(12) the Lorentz condition div(A) = takes the form 

d r {r 2 sin$g rr A r ) + d#{r 2 sin$g M A#) = 

which can be rewritten as 

^(sintfA?) + sintf<9 r (r 2 A r ) = , (13) 

wherefrom it follows cot ■dr p T + d r (r 2 A r ) = while the latter entails 

C cot -dr^T 



„2 



(14) 



with a constant matrix C. Then we can see that it should put C = T = or 
else A r will not be spherically symmetric and the confining one where only the 
powers of r equal to ±1 are admissible. As a result we come to the conclusion 
that one should put A r = A$ = in (12). After this we have ([, ] signifies 
matrix commutator) F = dA + gA A A = -pr^Tdt Adr + vr v ~ x Adr Adip + 
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gr» +u [r, A]dt A dp which entails [with the help of (10)] 



vr v ~ x A ar v+v 

*F = sin -&rd& A dp — dt A dd - - W, A]dr A d& 

sin V sin 17 

and the Yang-Mills equations (3) turn into 

fi{fi + I)?* sin 2 #r = g 2 r^ +2U [A, [r, A]] , 

v[v - l)r u ~ 2 A = g 2 r 2 ^ +u [r, [r, A]] . (15) 

It is now not complicated to enumerate possibilities for obtaining the confining 
solutions in accordance with (15), where we should put ji = —1, v — 1. 

(1) r=0 or A—0. This situation does obviously not correspond to a confining 
solution 

(2) r = C A with some constant C . This case conforms to that all the 
parameters a a describing electric colour Coulomb field (see Subsection 
2.3) and the ones f3 a for linear magnetic colour field are proportional - 
the situation is not quite clear from physical point of view 

(3) Matrices r, A are not equal to zero simultaneously and both matrices 
belong to Cartan subalgebra of SU(3)-Lie algebra. The parameters a a , (3 a 
of electric and magnetic colour fields are not connected and arbitrary, i. 
e. they should be chosen from experimental data. The given situation 
is the most adequate to the physics in question and the corresponding 
confining solution is in essence the same which has been obtained as far 
back as in Ref. [1] from other considerations. 

One can slightly generalize the starting ansatz (12) taking it in the form 
A = (r^T + A')dt + (r u A + B')d<p with matrices A' = A a \ a , B' = B a \ a and 
constants A a , B a . Then considerations along the same above lines draw the 
conclusion that the nontrivial confining solution is described by r, A, A', B' 
belonging to Cartan subalgebra. 

Let us remind that, by definition, a Cartan subalgebra is a maximal abelian 
subalgebra in the corresponding Lie algebra, i. e., the commutator for any 
two matrices of the Cartan subalgebra is equal to zero (see, e.g., Ref. [8]). 
For SU(3)-Lie algebra the conforming Cartan subalgebra is generated by the 
Gell-Mann matrices A3, A 8 which are 



A, 



(I 0\ , / 1 ON 



0-10 
\o 0; 



1 

\0 -2 J 



(16) 



Then it is easy to gain the only nontrivial solution in question in the form 
(which reflects the fact that for any matrix T from SU(3)-Lie algebra we have 
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Tr T = 0) 

Al+±=A\ = -^+A,, -Al+^Al = -^+A 2 , -4=A* = ^±^-(A 1+ A 2 ), 
W3 r V3 r r 

K+^K = b ir+B 1} -K+^K = b2r+B 2 , ~-j=A% = -(bi+b 2 )r-(B 1 +B 2 ), 

(17) 

where real constants dj, Aj, bj, Bj parametrize the solution, and we wrote down 
the solution in the combinations that are just needed to insert into the corre- 
sponding Dirac equation (see Section 3). From here it follows one more form 

Al = [(a 2 - ai)/r + A ± - A 2 ]/2, A s t = [A, + A 2 - ( Gl + a 2 )/r]VS/2 , 

K = [(bi ~ b 2)r + B,- B 2 }/2, Al = [(6i + b 2 )r + B 1 + B 2 }V3/2 (18) 

Clearly, the obtained results may be extended over all SU(iV)-groups with 
N > 2 and even over all semisimple compact Lie groups since for them the 
corresponding Lie algebras possess just the only Cartan subalgebra. Also we 
can talk about the compact non-semisimple groups, for example, U(iV). In 
the latter case additionally to Cartan subalgebra we have centrum consisting 
from the matrices of the form al^ {In is the unit matrix iV x N) with ar- 
bitrary constant a. The most relevant physical cases are of course U(l)- and 
SU(3)-ones (QED and QCD), therefore we shall not consider further general- 
izations of the results obtained but let us write out the corresponding solution 
for U(l)-case which will useful for to interpret above solutions (17)-(18) in 
the more habitual physical terms. It should also be noted that the nontrivial 
confining solutions obtained exist at any gauge coupling constant g, i. e. they 
are essentially nonperturbative ones. 



2.3 U(l)-case 

Under this situation the Yang-Mills equations (3) turn into the second pair of 
Maxwell equations 

d*F = (19) 
with F = dA, A = A^dx^ . We search for the solution of (19) in the form 
A = A t (r)dt + A^^dtp. It is then easy to check that F = dA = —d r A t dt A 
dr + dr-A^dr A dip and according to (10) we get *F = r 2 middrAtdd A dip — 
-^^drA^dt A d-d. From here it follows that (19) yields 

d r (r 2 d r A t ) = 0, d 2 r A v = , (20) 

and we write down the solutions of (20) as 

A t = - + A,A^ = br + B (21) 
r 
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with some constants a, b, A, B parametrizing solutions (further for the sake of 
simplicity let us put a = 1, b = 1 GeV, A = B = 0). 



To interpret solutions (21) in the more habitual physical terms let us pass on 
to Cartesian coordinates employing the relations 

dip dip 

ip = arctan(y/:r), dip = -—dx + -—dy (22) 

ox dy 

which entails 

Apdip = Ty dx + J X dy (23) 

x 2 + y l x 2 + y 2 

and we conclude that the solutions of (21) describe the combination of the 
electric Coulomb field with potential $ = A t and the constant magnetic field 
with vector-potential 

A = (A„ A s , A,) = (-^, 0) = (-^, ^, 0) , (24) 

x 2 + y 2 x 2 + y 2 sin v sin w 

which is linear in r in spherical coordinates and the 3-dimensional divergence 
divA = 0, as can be checked directly. Then Eqs. (19) in Cartesian coordinates 
takes the form 

A$ = 0, rotrotA = AA = (25) 

with the Laplace operator A = d 2 x + dy + d\. At last, it is easy to check that 
the solution under consideration satisfies the Lorentz condition div(A) = 0. 

Practically the same considerations as the above ones in electrodynamics show 
that the solutions (17)-(18) describe the configuration of the electric Coulomb- 
like colour field (components A t ) with potentials $ 3 , $ 8 and the constant mag- 
netic colour field (components A v ) with vector-potentials A 3 , A 8 which are lin- 
ear in r in spherical coordinates with 3-dimensional divergences divA 3 ' 8 = 
and A$ 3 ' 8 = 0, rotrotA 3 ' 8 = AA 3 ' 8 = 0. 



3 Relationship with QCD 



The previous section leads us to the following problem: how to describe possi- 
ble relativistic bound states in the obtained confining SU(3)-Yang-Mills fields? 
The sought description should be obviously based on the QCD-Lagrangian. 
Let us write down this Lagrangian (for one flavour) in arbitrary curvilinear 
(local) coordinates in Minkowski spacetime 

C = - fi o m - \{g^g^ - g^g^F^F^, n<v,a<(3 (26) 
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where, if denoting S(M) and £, respectively, the standard spinor bundle and 
3-dimensional vector one (equipped with a SU(3)-connection with the corre- 
sponding connection and curvature matrices A = A^dx^ = A^Xadx^, F = 
dA + gA A A = F^Xadx^ A dx v ) over Minkowski spacetime, we can con- 
struct tensorial product S = S(M) £g> £. It is clear that \& is just a section of 
the latter bundle, i. e. ^ can be chosen in the form ^ = (^1,^/2,^3) with 
the four-dimensional Dirac spinors fyj representing the j-th colour component 
while \& = ^(7° (g) I 3 ) is the adjont spinor, (f) stands for hermitian conjuga- 
tion, ® means tensorial product of matrices, (j, is a mass parameter, £> is the 
Dirac operator with coefficients in £ (see below). At last, we have the condition 
Tr(A a Afe) = 25 a b so that the third addendum in (26) has the form G(F,F)/8 
with G of (11), where coefficient 1/8 is chosen from physical considerations. 

From general considerations (see, e. g., Ref. [9]) the explicit form of the op- 
erator V in local coordinates x^ on Minkowski spacetime can be written as 
follows 

V = *( 7 e ® h)E» (dp ®h- ^^tV ®h~ igA^ ,a<b, (27) 

where the forms oj a b = oj^dx 11 obey the Cartan structure equations de a = 
uj a b Ae b , while the orthonormal basis e a = e a ^dx^ in cotangent bundle and dual 
basis E a = E^d^ in tangent bundle are connected by the relations e a {Ef ) ) = 
5%. At last, matrices r y a represent the Clifford algebra of the quadratic form 
Qi,z = x o — x 1 _ x \ ~ x \ i n C 4 . It should be noted that Greek indices //, u, ... are 
raised and lowered with g^ of (1) or its inverse g^ v and Latin indices a,b, ... 
are raised and lowered by r] a b = i] ab = diag(l,-l,-l,-l) except for Latin indices 
connected with SU(3)-Lie algebras, so that e^e b g^ u = rj ab , E^E^g^ = r\ ah and 
so on but A a = A a . 

Under the circumstances we can obtain the following equations according to 
the standard prescription of Lagrange approach from (26) 

W = hi* , (28) 

d* F = g(*F A A — A A *F) + gj , (29) 
where the source J (a nonabelian SU(3)-current) is 

j = 3 ;\ a * {dx») = *j = *u;\ a dx») = *Cf a„) (30) 

with currents 

f = j-dx" = #( 7/ , ® h)\ a * dx» , 
so summing over a = 1, 8 is implied in (26) and (30). 

When using the relation (see, e. g. Refs. [10]) r ) c E^<jJ liah ^ a ^ b = oJ fia i ) ^^^ a ^ b = 
— div(7) with matrix 1-form 7 = ^^dx^ [where div is defined by relation (4)] 
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and also the fact that (7^7° = 7°7 M , the Dirac equation for spinor \& will be 
id^(r <g> J 3 ) + ^div( 7 ) <g> h ~ g^(Y ® / 3 )^A a = -/x * • (28') 

Then multiplying (28) by ^A a from left and (28') by A a \l/ from right and 
adding the obtained equations, we get div(j a ) = div(j) = if spinor \& obeys 
the Dirac equation (28). 

The question now is how to connect the sought relativistic bound states with 
the system (28)-(29). To understand it let us apply to the experience related 
with QED. In the latter case Lagrangian looks like (26) with changing group 
SU(3)— >U(1) so ^ will be just a four-dimensional Dirac spinor. Then, as is 
known (see, e. g. Ref. [11]), when passing on to the nonrelativistic limit the 
Dirac equation (28) converts into the Pauli equation and further, if neglect- 
ing the particle spin, into the Schrodinger equation, parameter /x becoming 
the reduced mass of two-body system. The modulo square integrable solu- 
tions of the Schrodinger equation just describe bound states of a particle with 
mass fiQ, or, that is equivalent, of the corresponding two-body system. His- 
torically, however, everything was just vice versa. At first there appeared the 
Schrodinger equation, then the Pauli and Dirac ones and only then the QED 
Lagrangian. In its turn, possibility of writing two-body Schrodinger equation 
on the whole owed to the fact that the corresponding two-body problem in 
classical nonrelativistic (newtonian) mechanics was well posed and actually 
quantizing the latter gave two-body Schrodinger equation. Another matter 
was Dirac equation. Up to now nobody can say what two-body problem in 
classical relativistic (einsteinian) mechanics could correspond to Dirac equa- 
tion. The fact is that the two-body problem in classical relativistic mechanics 
has so far no single-valued statement. Conventionally, therefore, Dirac equa- 
tion in QED is treated as the relativistic wave equation describing one particle 
with spin one half in an external electromagnetic field. 

There is, however, one important exclusion - the hydrogen atom. When solv- 
ing the Dirac equation here one considers mass parameter /xo to be equal to 
the electron mass and one gets the so-called Sommerfeld formula for hydrogen 
atom levels which passes on to the standard Schrodinger formula for hydrogen 
atom spectrum in nonrelativistic limit (for more details see, e. g. Ref. [11]). 
But in the Schrodinger formula mass parameter /xq is equal to the reduced 
mass of electron and proton. As a consequence, it is tacitly supposed that in 
Dirac equation the mass parameter should be equal to the same reduced mass 
of electron and proton as in Schrodinger equation. Just the mentioned reduced 
mass is approximately equal to that of electron but, exactly speaking, it is not 
the case. We remind that for the problem under discussion (hydrogen atom) 
the external field is the Coulomb electric one between electron and proton, 
essentially nonrelativistic object in the sense that it does not vanish in non- 
relativistic limit at c — > 00. If now to place hydrogen atom in a magnetic field 
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then obviously spectrum of bound states will also depend on parameters de- 
cribing the magnetic field. The latter, however, is essentially relativistic object 
and vanishes at c — > oo because, as is well known, in the world with c = oo 
there exist no magnetic fields (see any elementary textbook on physics, e. g. 
Ref. [12]). But it is clear that spectrum should as before depend on fx as well 
and we can see that /io is the same reduced mass as before since in nonrela- 
tivistic limit we again should come to the hydrogen atom spectrum with the 
reduced mass. So we can draw the conclusion that if an electromagnetic field 
is a combination of electric Coulomb field between two charged elementary 
particles and some magnetic field (which may be generated by the particles 
themselves) then there are certain grounds to consider the given (quantum) 
two-body problem to be equivalent to the one of motion for one particle with 
usual reduced mass in the mentioned electromagnetic field. As a result, we can 
use the Dirac equation for finding possible relativistic bound states for such 
a particle implying that this is really some description of the corresponding 
two-body problem. 

Actually in QED the situaion is just as the described one but magnetic field 
is usually weak and one may restrict oneself to some corrections from this 
field to the nonrelativistic Coulomb spectrum (e. g., in the Seemann effect). 
If the magnetic field is strong then one should solve just Dirac equation in 
a nonperturbative way (see, e.g. Ref. [13]). The latter situation seems to be 
natural in QCD where the corresponding (colour) magnetic field should be very 
strong (see Section 4) because just it provides linear confinement of quarks as 
we could see above (see also Refs. [2-4]) and the given field also vanishes in 
nonrelativistic limit (for more details see Refs. [2,3]). 

At last, we should make an important point that in QED the mentioned elec- 
tromagnetic field is by definition always a solution of the Maxwell equations 
so within QCD we should require the confining SU(3)-field to be a solution 
of Yang-Mills equations. Consequently, returning to the system (28)-(29), we 
can suggest to decribe relativistic bound states of two quarks (mesons) in 
QCD by the compatible solutions of the given system. To be more precise, 
the meson wave functions should be the nonperturbative modulo square inte- 
grable solutions of Dirac equation (28) (with the above reduced mass /i ) in 
the confining SU(3)-Yang-Mills field being a nonperturbative solution of (29). 
In general case, however, the analysis of (29) is difficult because of availabil- 
ity of the nonabelian current J of (30) in the right-hand side of (29) but we 
may use the circumstance that the corresponding modulo square integrable 
solutions of Dirac equation (28) might consist from the components of form 
tyj ~ r Q ^e~ /3jr with some aj > 0, f3j > which entails all the components of 
the current J to be modulo << 1 at each point of Minkowski space. The latter 
will allow us to put J « and we shall come to the problem of finding the con- 
fining solutions for the Yang-Mills equations of (29) with J = whose unique 
nontrivial form has been described in previous Section and after inserting the 
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found solutions into Dirac equation (28) we should require the corresponding 
solutions of Dirac equation to have the above necessary behaviour. Under the 
circumstances the problem becomes self-consistent and can be analyzable and 
this has been done actually in Ref. [1]. 

It is clear that all the above considerations can be justified only by comparison 
with experimental data but now we obtain some intelligible programme of 
further activity which has been partly realized in Refs. [2-4]. 



4 Verification of confinement criterion and estimates of colour mag- 
netic field strength 

To illustrate some of the above let us verify the Wilson confinement criterion 
[14] for the confining solutions (17)— (18) and also adduce numerical estimates 
for strength of colour magnetic field responsible for linear confinement in the 
ground state of charmonium. 



4-1 Verification of confinement criterion 

The Wilson confinement criterion is in essence the assertion that the so-called 
Wilson loop W(c) should be subject to the area law for the confining glu- 
onic field configuration. In its turn, the latter law is equivalent to the fact 
that energy E(R) of the mentioned configuration (gluon condensate) is lin- 
early increasing with R, a characteristic size of some volume V containing 
the condensate. We can easily evaluate E{R) for solutions (17)— (18) using the 
Too-component (volumetric energy density) of the energy-momentum tensor 



It is not complicated to obtain the curvature matrix (field strentgh) corre- 
sponding to the solution (17)-(18) 



F = F a Xadx^ A dx v = -d r (A a t \ a )dt Adr + d r (A a \ a )dr A dip , (32) 



for a SU(3)-Yang-Mills field 



T = — ( — 

47T I 



pa pa a/3 , lp pa . a/3 ' j8 



) 



(31) 



which entails the only nonzero components 




di — CL2 




(ai + a 2 )VS 




b\ ~ h 
2 




2 



(33) 



2r 2 



2r 2 
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and, in its turn, 



Too = In 



1 

An 



trJ 



+ 



1 



4r 2 sin 2 d 



(Kir + (K) 



8 \2 



rip 



Tip) 



a 



a\a 2 + a\ b\ 

A ~t~ 



b x b 2 + b\ 



A 



B 



16n \ r 4 3r 2 sin 2 $ J r 4 ^ r 2 sin 2 -(9 ^ ^ 

with „4 > 0, B > 0. Let V be the volume between two concentric spheres 
with radii R < R. Then E(R) = J v T 00 VSd 3 x = J v T 00 r 2 sin ddrdfidip . As 
was discussed in Ref. [4], the notion of classical gluonic field is applicable only 
at distances » Xb, the de Broglie wavelength of quark, so we may take Rq 
to be of order of a characteristic size of hadron (meson). Besides it should 
be noted that classical T 00 of (34) has a singularity along ,2-axis (■# = 0, n) 
and we have to introduce some angle -#o whose physical meaning is to be 
clarified a little below. Under the circumstances, with employing the relations 
/ dd/ sin'i? = lntan-^/2, tan^/2 = sin-^/ (1 + cos-#) = (1 — cosi?)/ sin-^, we 
shall have 

- / 77(0 + d*) -™ - 

(35) 

with E = ^-2ttBRo In It is clear that at i? -> oo we get E(R) ~ 

i. e. the area law is fulfilled. 



4-2 Estimates of colour magnetic field strength 



To estimate R and $0 under real physical situation let us consider the ground 
state of charmonium rj c (lS) and apply the ideology of Section 3 for describ- 
ing the given state. We use the parametrization of relativistic spectrum of 
charmonium from Ref. [4]. Namely, Table 1 contains necessary parameters of 
solutions (17)— (18) for the sought estimates. As for parameters A 12 of (17)- 

Table 1 



Gauge coupling constant, mass parameter no and parameters of the confining SU(3)- 
gluonic field for charmonium. 



9 


Mo 
(GeV) 


a\ 


a 2 


bi 
(GeV) 


b 2 
(GeV) 




B 2 


0.46900 


0.62500 


2.21104 


-0.751317 


20.2395 


-12.6317 


6.89659 


6.89659 



(18) then we can put them equal to zero [4]. One can note that the reduced 
mass parameter fi Q is consistent with the present-day experimental limits [15] 
where the current mass of c-quark (2/i ) is accepted between 1.1 GeV and 1.4 
GeV. Then wave function \P = (vE^, \J/ 2? ^3) of the state under consideration 
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[as a modulo square integrable solution of Dirac equation (28) in the field 
(17)— (18)] is given by the form (with Pauli matrix a±) 

•^^■-■(^X))^ 1 - 2 - 8 (36) 

with the 2D eigenspinor $j = ^ ^ of the euclidean Dirac operator on the 

unit sphere § 2 . The explicit form of $j is not needed here and can be found 
in Ref. [4]. The radial parts are 



P. 




with aj = yJA 2 - g 2 a 2 , Aj = -1 - gBj, (5j = yj $ - u] + g 2 b 2 , Pj = gbj + fy, 
Qj = + Vj where a 3 = —(a± + a 2 ),b 3 = — (&i + b 2 ),B 3 = —(B 1 + B 2 ) and 
constants Cj are determined from the normalization conditions Jo°°(l-^ji| 2 + 
\Fj 2 \ 2 )dr = 1/3. Energy of the ground state is obtained as 

JU _ -Ai5f 2 ai&i + ai|Ai|/i , -A 2 g 2 a 2 b 2 + a 2 \A 2 \no -A 3 g 2 a 3 b 3 - a 3 \A 3 \/i n _ n „ - r 
e - Z^i = 72 1 72 1 72 - ^ y ' yb UeV • 

j=l A l A 2 A 3 

(38) 

Then we can put R = |E; 3 =il/^ « 0.0399766 fm and, adding the rest 
energies of quarks 2(2/z ) to E of (35), we consider e = -E +4/i = 2.9796 GeV 
which entails i?o ~ 0.723388 and one may consider $ to be the parameter 
determining some cone < i? < $o so the quark emits gluons outside of the 
cone. It is evident that the Wilson criterion is purely classic one: quantum 
mechanically the stay probability of quarks at distances > Ro is of order 
Y, 3 j=i(\Fji\ 2 + |-Fj 2 | 2 ) and is damped exponentially fast at R > R . Just the 
colour magnetic field defines this damping through the coefficients f3j. 

At last, using the Hodge star operator in 3-dimensional euclidean space [where 
ds 2 = g^dx* <g) dx v = dr 2 + r 2 (dd 2 + sin 2 tfdtp 2 ), *(dr A dtff) = sin tidtp, *(dr A 
dtp) = — cM/sini?, *(d$ A dtp) = drjir 2 sin •&], we can confront components 
F^, F^ v of (33) with 3-dimensional 1-forms of the magnetic colour field 

2 sin t? 2 sin 



I 38 38 

that are modulo equal to y g^H^ Hv or 



R3 = \h-hl ff8 = Ih + hlVl 

2rsin$' 2rsin$ ' v ; 

where we shall put sin-$ = 1 for simplicity. Table 2 contains estimates of 
magnetic colour field strength in the ground state of charmonium whereas the 
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Bohr radius ao = 0.529177249 • 10 5 fm [15]. Also when calculating we applied 
the relations 1 GeV" 1 « 0.21030893 f m , 1 T « 0.692508 x 10~ 15 GeV 2 . 

Table 2 

Magnetic colour field strengths in the ground state of charmonium. 



n CI <?V 


Ro = 0.0399766 fm 






r 






H 8 




(fm) 


(T) 




(T) 




O.IRo 


0.124857 


10 19 


0.500516 


10 18 


Rq 


0.124857 


10 18 


0.500516 


10 17 


WR 


0.124857 


10 17 


0.500516 


10 16 


1.0 


0.499137 


10 16 


0.200089 


10 16 


ao 


0.943232 


10 11 


0.378114 


10 11 



It is seen that at the characteristic scales of the charmonium ground state the 
strength of magnetic colour field responsible for linear confinement reaches 
huge values of order 10 17 -10 18 T. For comparison one should notice that the 
most strong magnetic fields known at present have been discovered in magnetic 
neutron stars, pulsars (see, e. g., Ref. [16]) where the corresponding strengths 
can be of order 10 9 -10 10 T. So the characterictic feature of confinement is 
really very strong magnetic colour field between quarks. In a certain sense the 
essence of confinement can be said to be just in enormous gluon concentrations 
and magnetic colour field strentghs in space around quarks. 



5 Concluding remarks 

The aim of this Letter we pursued was to specify a scenario of linear confine- 
ment of quarks, at any rate, for mesons and quarkonia. As we could see, crucial 
step consisted in studying exact solutions of the SU(3)-Yang-Mills equations 
modelling confinement. Namely structure of the confining solutions for SU(3)- 
Yang-Mills equations prompts to the idea that linear confinement should be 
governed by magnetic colour field. But the latter is essentially relativistic ob- 
ject vanishing in nonrelativistic limit which implies the relativistic effects to 
be extremely important for the confinement mechanism [2-4]. 

Further specification of the confinement mechanism can consist in the in- 
quiry answer: what the way of gluon exchange between quarks is? The matter 
concerns such a modification of gluon propagator which might correspond to 
linear confinement at large distances. As was remarked in Ref. [2] , the sought 
modification should be carried out on the basis of the exact solutions of the 
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SU(3)-Yang-Mills equations modelling quark confinement. We hope to con- 
sider this modification elsewhere. 

There is also an interesting possibility of indirect experimental verification of 
the confinement mechanism under discussion. Really solutions (21) point out 
the confinement phase could be in electrodynamics as well. Though there exist 
no elementary charged particles generating a constant magnetic field linear in 
r, the distance from particle, after all, if it could generate this elecromagnetic 
field configuration in laboratory then one might study motion of macroscopic 
charged particles in that field. The confining properties of the mentioned field 
should be displayed at classical level too but the exact behaviour of particles 
in this field requires certain analysis of the corresponding classical equations 
of motion. 
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